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PREFACE 

Laser  Doppler  velocimetry  has  been  under  development  for  approxi- 
mately ten  years.  Initially,  laser  Doppler  velocimetry  was  developed  as  a 
research  tool.  On  25  June  1975,  Eastern  Airline  Flight  66  crashed  while 
making  an  approach  to  John  F.  Kennedy  International  Airport  in  New  York. 

The  cause  of  the  crash  has  been  attributed  to  wind  shear,  and  a significant 
program  in  wind  shear  detection  resulted.  The  laser  Doppler  velocimeter  is 
a candidate  for  the  remote  measurement  of  atmospheric  wind  and,  thus,  the 
detection  of  wind  shear.  In  order  to  assess  the  ability  of  the  laser  Doppler 
velocimeter  to  measure  winds  remotely,  a test  for  the  comparison  of  wind 
measured  by  the  laser  Doppler  velocimeter  and  wind  measured  by  tower- 
mounted  anemometers  was  conducted  at  the  National  Oceanic  and  Atmospheric 
Administration's  Table  Mountain  test  site  near  Boulder,  Colorado.  The  laser- 
measured  winds  showed  good  agreement  with  the  anemometer -measured  winds. 
However,  because  of  the  limited  altitude  of  the  tower,  verification  was  limited 
to  altitudes  less  than  150  meters  In  order  to  validate  the  laser  Doppler 
velocimeter  at  altitudes  above  150  meters,  a test  similar  to  the  Table  Moun- 
tain test  was  conducted  near  the  481 -meter  instrumented  WKY-TV  television 
transmission  tower  at  the  National  Severe  Storms  Laboratory  test  site  near 
Norman,  Oklahoma  This  report  describes  the  results  of  that  test. 


The  assistance  of  R.  Craig  Goff  of  the  National  Severe  Storms  Laboratory 
in  the  performance  of  the  tests  is  gratefully  acknowledged. 
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1.  INTRODUCTION 


I.)  BACKGROUND 

Significant  effort  is  currently  being  devoted  to  the  development  of  in- 
strumentation to  remotely  sense  atmospheric  flow  phenomena.  Some  of  the 
avenues  being  pursued  are  active  and  passive  acoustic  sensors,  optical 
sensors,  and  radio  methods.  A useful  survey  of  such  methods  is  presented 
in  Ref  1.  Two  advantages  of  remote  sensors  are  that  flow  conditions  can 
be  ascertained  in  regions  of  space  where  it  would  not  be  convenient  to  locate 
instrumentation  hardware,  and  no  interference  with  the  flow  at  the  point  of 
interest  is  introduced  by  their  use.  The  laser  Doppler  velocimeter  (LDV) 
is  a particularly  attractive  device  for  remote  sensing  of  atmospheric  phe- 
nomena. In  the  LDV  system,  the  laser  radiation  backscattered  by  moving 
particulates  in  the  atmosphere  is  used  to  determine  the  velocity  of  the  flow. 
Since  it  is  possible  to  direct  the  laser  focal  volume  at  a selected  sequence 
of  points  in  space,  data  from  a scanning  LDV  system  can  be  used  to  deter- 
mine the  velocity  field  rapidly  and  over  a range  of  altitudes.  A CO^  laser 
Doppler  velocimeter  system  has  the  following  advantages  over  other  remote 
sensing  techniques:  (1)  the  sensing  volume  can  be  varied  with  ease  as  only 
optic  pointing  and  focusing  operations  are  involved;  (2)  the  ambient  aerosol 
provides  a sufficient  scattering  target;  and  (3)  the  sensing  mechanism  is  non- 
mechanical, which  results  in  the  potential  for  a high  frequency  turbulence 
sensor. 

The  feasibility  of  using  an  LDV  system  for  the  remote  sensing  of 
low  altitude  winds  and  for  the  detection  and  tracking  of  aircraft  wake  vortices 
has  been  demonstrated  (Refs.  2,  3, and  4).  However,  the  development  of  an 
effective  LDV  system  for  monitoring  wind,  wind  shear,  and  wake  vortices 
required  a further  refinement  and  application  of  the  technology. 
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Previous  field  testing  of  the  LDV  system  was  accomplished  at  the 
John  F.  Kennedy  International  Airport  (JFK)  and  at  the  Wave  Propagation 
Laboratory  of  the  National  Oceanic  and  Atmospheric  Administration  near 
Boulder,  Colorado  (Refs.  5 and  6).  Because  of  the  demonstrated  capability 
to  measure  winds  at  these  test  sites,  it  was  deemed  appropriate  to  assess 
the  ability  of  the  laser  Doppler  velocimeter  to  measure  winds  at  higher 
altitudes. 

1.2  PROGRAM  OBJECTIVES 

Thi  primary  objective  of  the  tests  at  the  National  Severe  Storms  Lab- 
oratory was  to  verify  the  ability  of  the  laser  Doppler  velocimeter  to  remotely 
measure  winds  at  altitudes  up  to  450  m.  Additionally,  it  was  desired  to  mea- 
sure wind  phenomena  in  thunderstorms  which  were  expected  to  pass  over  the 
test  area  during  the  test  period.  The  second  objective  was  only  partially  ac- 
complished as  only  one  thunderstorm  passed  near  the  test  area.  Its  intensity 
was  not  sufficient  to  assess  the  measurement  of  wind  phenomena  in  severe 
storms.  However,  data  were  obtained  in  rain,  and  the  results  are  presented 
herein. 

1.3  REPORT  FORMAT 

A brief  description  of  the  mobile  laser  Doppler  system  is  presented 
in  the  following  sections.  A more  complete  description  is  presented  in 
Ref.  5.  A description  of  the  LDV  contained  in  the  Lockheed  Mobile  Atmos- 
pheric Unit  (MAU)  is  presented  in  Section  2 followed  by  a description  of  the 
computer  software  algorithms  in  Section  3 The  field  tests  at  the  National 
Severe  Storms  Laboratory  site  are  described  in  Section  4,  along  with  analysis 
of  the  LDV  data  and  comparisons  with  the  meteorological  tower  data.  Section 
5 presents  the  conclusions  and  recommendations. 


2.  LASER  DOPPLER  SYSTEM  DEVELOPMENT 


2.1  SYSTEM  DESCRIPTION 

An  LDV  remote  wind  sensor  senses  air  movement  by  measurement  of 
the  Doppler  frequency  shift  of  laser  radiation  backscattered  by  the  atmos- 
pheric aerosol.  An  instrument  must  incorporate  means  to  transmit  the  laser 
radiation  to  the  region  of  interest,  collect  the  radiation  scattered  from  the 
atmospheric  aerosol,  and  to  photomix  on  a photodetector  the  scattered  radia- 
tion and  a portion  of  the  transmitted  beam.  The  difference  between  the  trans- 
mitted frequency  and  the  returned  frequency  is  the  Doppler  shift  frequency. 

The  Doppler  frequency  shift  signal  is  generated  at  the  photodetector  and  is 
directly  proportional  to  the  magnitude  of  the  wind-velocity  component  in  the 
direction  of  the  line-of-sight  of  the  laser  beam.  This  velocity  component  is 
hereafter  called  the  line-of-sight  velocity.  The  magnitude  of  the  Doppler  shift, 
Af,  is  given  by 

Af  = | v|  cos9, 

where 

V = the  velocity  vector  in  the  region  being  sensed, 

6 = the  angle  subtended  by  the  velocity  vector  and 
the  optic  system  line-of-sight,  and 

X = the  laser  radiation  wavelength  (10.6  microns 
for  CO^  laser). 

A Doppler  shift  of  188  kHz  results  per  m/sec  of  line-of-sight  velocity  com- 
ponent. Measurement  of  the  Doppler  shift  frequency,  Af,  yields  directly  the  line- 
of-sight  velocity  component,! V|  cos0.  Some  typical  advantages  of  the  laser 
Doppler  method  are:  (1)  the  Doppler  shift  is  a direct  absolute  measure  of  the 
line-of-sight  velocity  (for  example,  a hot  wire  yields  wind  speed  via  a cooling 
effect  on  the  wire);  (2)  the  ease  with  which  the  sensing  volume  can  be  varied 
(only  optics  pointing  and  focusing  operations  being  involved);  (3)  the  ambient 

^Techniques  for  resolving  the  sign  of  the  line-of-sight  velocity  component  are 
available.  They  are  not  discussed  here  because  they  were  not  available  for 
the  subject  test. 
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aerosol  provides  sufficient  scattering,  thus  enabling  operation  in  "clear  air" 
conditions;  and  (4)  the  ambient  aerosol  tracer  has  a small  inertia  and  re- 
sponds quickly  to  variations  in  airspeed  and,  thus,  can  be  a good  turbulence 
indicator . 


A useful  instrument  must  also  incorporate  means  for  scanning  the 
focal  volume  through  a selected  sequence  of  points  in  space  and  to  effect 
the  required  signal  processing,  on-line  readout,  and  permanent  recording 
requirements.  The  hardware  implementation  of  the  mobile  laser  Doppler 
unit  used  during  this  investigation  is  discussed  in  the  following  sub- 
sections. The  overall  configuration  is  summarized  in  Fig.  2-1. 


2.1.1  Basic  Laser  Doppler  Velocimeter  Optical  System 

The  optical  system  is  of  monostatic  design  and  utilizes  a continuous- 
wave  las er . The  arrangement  depends  on  focusing  the  transmitter  telescope 
at  the  location  of  interest  for  its  spatial  resolution  property.  Details  of  the 
optical  arrangement  are  shown  in  Fig.  2-2. 


A horizontally  polarized,  20 -watt,  continuous-wave  CO^  laser  beam 
(10.6-micron  wavelength)  emerges  from  the  laser  (T)  and  is  deflected  90 
degrees  by  a mirror  (T)  and  by  a 90%  reflecting  beamsplitter  (7).  The  ap- 
proximately 6-mm  diameter  beam  then  passes  through  a Brewster  window  (7) 
and  a CdS  quarter  waveplate  (T)  which  converts  the  beam  to  circular  polari- 
zation. The  beam  impinges  on  the  secondary  mirror  (T)and  is  expanded  and 
reflected  into  the  primary  mirror  (30-cm  diameter)  (7)and  then  focused  oui 
into  the  atmosphere.  A small  portion  of  the  original  laser  beam  is  transmitted 
through  the  beamsplitter  (7)  and  is  used  as  a local  oscillator  after  being  ro- 
tated to  vertical  polarization  by  a half  waveplate  (9).  Energy  scattered  by 
aerosols  at  the  focal  volume  (T5)  is  collected  by  the  primary  mirror  (7),  col- 
limated by  the  secondary  (T),  and  passed  through  the  quarter  waveplate  (^. 

A wire  stop  (l6)  eliminates  most  of  the  secondary  mirror  reflection  of  the  out- 
going beam.  The  quarter  waveplate  changes  the  polarization  of  the  aerosol  back- 
scattered  radiation  from  circular  to  vertical  linear  polarization.  The  vertically 
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FIGURE  2-2.  TYPICAL  OPTICAL  COMPONENT  CONFIGURATION  OF  LOCKHEED 
LASER  DOPPLER  VELOCIMETER. 


polarized  beam  is  approximately  78%  reflected  off  the  Brewster  window  (?)  and 
directed  through  the  beamsplitter  (o)  where  it  is  combined  with  the  local  oscil- 
lator radiation.  After  passing  through  the  collecting  lens  (fl) , the  two  beams 
are  photomixed  on  the  detector  (?2)  in  a heterodyne  configuration.  The  electrical 
output  of  the  detector  (£2)  is  amplified  (l?)  with  a 5-MHz  bandwidth,  20-dB  gain 
low -noise-type  preamplifier  and  fed  into  a spectrum  analyzer  (Ti)  . 

An  alternative  operating  configuration  consisted  of  utilizing  the  portion 
of  the  outgoing  beam  backscattered  into  the  interferometer  by  the  secondary 
mirror  (?)  as  the  local  oscillator  beam.  This  mode  of  operation  is  less  sus- 
ceptible to  optical  misalignment  difficulties  and  was  the  technique  used 
during  this  investigation.  When  incorporated,  the  optical  leg  (T)(z)  (9)  was 
deactivated  by  the  beam  stop  (?)  and  the  wire  stop  (l6)  was  removed. 

2.1.2  Optical  Scanning  System 

In  order  to  provide  the  flexibility  to  operate  in  various  modes  for  which 
the  system  was  designed,  a scanning  arrangement  as  shown  in  Fig.  2-3  is  used. 
Modes  of  operation  include  vortex  tracking  (not  required  for  the  measurements 
described  in  this  document)  and  Velocity  Azimuth  Display  (VAD)  for  the  meas- 
urement of  atmospheric  wind.  The  mirror  assembly,  AB,  can  be  rotated  about 
the  vertical  axis  for  the  scanning  in  azimuth  necessary  for  the  VAD  (also  called 
the  conical-scan  mode  of  operation).  Mirror  A is  adjusted  to  control  the 
elevation  angle  of  the  beam,  thus  controlling  the  cone  angle  of  the  conical 
scan.  The  scanning  hardware  as  deployed  on  the  mobile  van  is  shown  in 
Fig. 2-4. 


Range  scanning  of  the  system's  focal  volume  is  accomplished  by  varying 
the  distance  between  the  telescope  secondary  mirror,  E,  and  the  primary  mir- 
ror, D.  This  is  effected  by  varying  the  position  of  the  mirror,  E,  in  a con- 
trolled manner  by  an  electric-motor/ optical-encoder  combination. 

The  operator  inputs  for  the  scanning  system  are  made  through  a con- 
trol panel  incorporating  thumbswitch  controls  and  light-emitting -diode 
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monitors.  The  system's  scan  capabilities  are  summarized  in  Fig.  2-5.  A 
view  of  tlie  mobile  unit  is  presented  in  Fig.  2-6.  Interior  views  of  the  mobile 
unit  are  shown  in  Figs.  2-7  and  2-8. 

2.1.3  Signal  Processing  System 

The  Doppler  frequency  shift  of  the  photodetector  output  is  determined 
through  the  use  of  a sampled  spectrum  analyzer  which  provides  frequency 
spectra  (intensity  of  returned  signal  as  a function  of  Dopplti  shift)  at  a rate 
of  70  signatures  per  second.  A typical  Doppler  wind  signature  is  shown  in 
Fig.  2-9.  To  yield  a line-of-sight  velocity  estimate,  a voltage  which  has  the 
same  time  behavior  as  the  Doppler  shift,  f^,  given  by  the  peak  of  the  spectrum 
is  generated. 

The  implementation  of  this  technique  is,  in  essence,  a recursive  com- 
parison method.  The  spectrum  analyzer  scan  is  driven  by  a sawtooth  voltage 
derived  from  a D/A  converter,  the  input  to  which  is  counterc locked  at  a con- 
stant rate,  hence  the  digital  output  of  the  counter  represents  frequency  on  a 
linear  scale.  At  each  new  count,  the  spectrum  a;  alyzer  output  is  converted 
to  a digital  representation  by  an  A/D  converter,  and  the  binary  number  repre- 
senting the  current  sample  is  compared  with  that  obtained  on  the  previous 
count.  If  the  current  sample  is  the  larger  of  the  two,  it  is  saved  by  storing 
in  a latch,  along  with  the  binary  number  representing  its  frequency;  if  it  is 
smaller,  the  previous  one  is  retained  until  the  next  comparison.  This  process 
is  continued  for  the  entire  sweep.  It  is  evident  that  the  number  remaining  in 
the  frequency -store  latch,  when  the  sweep  is  completed,  corresponds  to  the 
highest  signal  power  observed;  i.e.,  the  peak  of  the  spectrum.  At  the  end  of 
each  sweep,  the  new  peak  frequency  replaces  that  obtained  on  the  previous 
sweep.  An  example  of  the  output  is  shown  in  Fig.  2- 10  for  the  case  of  an  FM 
signal  of  center  frequency  2.0  MHz  (fQ)  modulated  sinusoidally  to  +200  kHz 
about  f at  a 5-Hz  rate.  The  raw  spectral  information  (output  of  the  spectrum 
analyzer)  is  also  made  available  to  the  Systems  Engineering  Laboratories  (SFI.) 
810A  data-logging  minicomputer  which  is  programmed  to  generate  its  own  esti- 
mate of  the  spectral  peak. 
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FIGURE  Z-o.  LOCKHEED -HUNTSVILLE  MOBILE  ATMOSPHERIC  UNIT. 


FIGURE  2-7.  INTERIOR  VIEW  OF  MOBILE  ATMOSPHERIC  UNIT  LOOKING  FORWARD. 
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FIGURE  2-8.  INTERIOR  VIEW  OF  MOBILE  ATMOSPHERIC  UNIT. 


(10.6  m/ sec) 


Frequency  (velocity) 


Frequency  shift  at  full  scale 
is  selected  by  operator 


FIGURE  2-9.  TYPICAL  LASER  DOPPLER  VELOCIMETER  WIND 
SIGNATURE  AS  DISPLAYED  BY  SPECTRUM  ANALYZER. 
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FIGURE  2-10.  OUTPUT  OF  SIGNAL  PROCESSOR  FOR  FRF.Ql  FNC Y - 
MODULATED  INPUT.'- 

* 

A provision  is  included  for  tracking  s ingle-s ideband  suppressed-carrier 
signals,  with  an  identification  of  upper  or  lower  s(deband  such  that  when  used 
in  conjunction  with  an  acousto-optic  modulator  the  unit  can  discriminate  the 
sign  of  the  Doppler  shift.  The  signal  feedthrough  at  the  translated  frequency 
can  also  be  discriminated  against  digitally,  thus  eliminating  the  need  for  a 
"notch  filter."  Such  a translating  device  was  not  used  in  the  data  collection 
described  herein. 

2.1.4  Data  Recording  and  Display 

The  primary  datS-ga^her ing  function  is  performed  by  the  SEL  810A 
general-purpose  minicomputer.  Data  acquired  by  the  Mobile  Atmospheric 
Unit  are  formatted  by  the  computer  software  and  stored  on  magnetic  tape  for 
subsequent  off-line  processing.  The  SEL  7-track  tape  control  and  magnetic 
tape  units  allow  digital  recording  of  data  at  800  bpi  at  45  ips.  The  data 
logged  by  the  computer  includes: 
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All  scan  volume  location  parameters 
"Mode  of  operation"  identifier 

The  instantaneous  line-of- sight  velocity  information 
The  Doppler  spectrum  peak  strength 
Full  spectrum  intensity  and  frequency  information 
A data-quality  identifier. 

Properties  of  the  Doppler  spectrum;  namely,  the  amplitude  and  frequency 
corresponding  to  the  spectral  peak,  are  obtained  as  a result  of  on-line  com- 
puter processing  exc  ept  for  the  frequency,  which  is  also  obtained  by  the 
spectral  peak  locator  (velocity  processor)  discussed  previously.  The  latter 
allows  some  flexibility  for  on-line  operator  displays  (see  below).  The  veloc- 
ity processor  estimate  of  the  instantaneous  line-of-sight  velocity,  updated  at 
a 70-Hz  rate,  is  available  in  analog  format  which  can  be  recorded  directly  on 
a stripchart  recorder,  an  option  which  is  extremely  useful  during  the  VAD 
mode  of  operation  for  monitoring  the  characte  ristic  profile. 

A view  of  the  computer  and  associated  LDV  electronics  is  shown  in 
Fig. 2-11. 

2.2  WINDS  ALOFT  SENSING 

Using  the  basic  system  outlined  previously  it  is  possible,  by  scanning 
operations,  to  determine  the  three-component  wind  field  at  any  altitude  be- 
tween 16  and  865  meters.  The  scanning  method  employed  is  commonly  re- 
ferred to  as  the  Velocity  Azimuth  Display  (VAD)  technique  and  was  first  used 
by  Lhermitte  and  Atlas  in  conjunction  with  a microwave  radar  (Ref.  7). 

The  telescope  is  focused  at  the  altitude  of  interest,  the  beam  being  di- 
rected at  a zenith  angle,  a.  The  beam  is  then  scanned  in  azimuth,  thus  tracing 
out  a circle  at  the  selected  altitude  (Fig.  2- 12). 
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FIGURE  Z-  1 1 . COMPUTER  MAINFRAME  TELETYPE  AND  LASER 
DOPPLER  VELOCIMETER  ELECTRONICS. 


Volume 
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FIGURE  2-12.  PRINCIPLE  OF  VELOCITY  AZIMUTH  DISPLAY  OPERATION. 


The  instantaneous  line-of-sight  component  of  velocity  within  the  sensing 
volume  as  measured  by  the  LDV,  vy,  is  given  by 

v = v,  sina  cos(0  - 0 ) + w cosa  , 
r h o 

v and  9 , respectively,  being  the  speed  and  direction  of  the  horizontal  wind 

h o’ 

component,  and  w the  vertical  component  at  the  height  being  sampled.  The 
azimuthal  dependence  of  v is  sufficient  to  yield  the  horizontal  speed  and 
direction  and  the  vertical  component  of  velocity. 

In  the  present  mode  of  operation,  the  system  is  unable  to  distinguish 
between  positive  and  negative  values  of  vr>  Therefore,  it  is  the  absolute 
value  of  vr  (|vr|)  that  is  sensed.  This  results  in  a signal  as  shown  in  Fig. 
2-13  instead  of  the  sinusoidal  signal  as  shown  in  Fig.  2-12.  Thus,  there  is 
an  ambiguity  of  180  deg  in  the  wind  direction  since  it  is  uncertain  which  peak 
in  Fig.  2- 13  represents  looking  into  the  wind.  In  practice,  the  ambiguity  is 
removed  by  the  operator  recording  approximate  wind  direction.  The  data 
processing  technique  can  then  calculate  exact  wind  direction.  This  resolves 
all  wind  direction  ambiguities  if  the  operator's  input  estimate  is  within  +89 
deg  of  the  true  wind  direction. 


Vh  = (V1  + V2)/,(2  Slna> 
w = (Vj  - V 2>/(2  cosa) 


FIGURE  2-13.  AZIMUTH  ANGLE  DEPENDENCE  OF  MEASURED 
VELOCITY  COMPONENT 
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While  operating  in  the  VAD  mode  the  system  is  capable  of  measuring 
winds  at  n (n  - 1 through  8)  altitudes  (th  t can  be  dialed  in  by  using  thumb- 
switches)  in  sequence  over  a total  time  period  of  5np  seconds  where 

5 sec  = time  for  conical  sweep  at  one  altitude  for  this  test, 

n = number  of  altitudes  to  be  interrogated,  and 

p - number  of  VAD  scans  at  each  altitude  (can  be  chosen 
to  be  1 through  7). 

During  this  investigation,  all  data  were  taken  in  the  single  altitude  mode  (i.e., 
n 1,  p = 1),  although  data  were  measured  at  several  different  altitudes  for  the 
various  recorded  runs  of  the  test. 
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3.  COMPUTER  SOFTWARE  SYSTEM  DEVELOPMENT 

Acquisition  and  processing  of  the  LDV  signature  are  accomplished  by 
means  of  a compact  data-handling  system  developed  specifically  for  the  Lock- 
heed"Hunts  ville  MAU.  The  general  elements  of  the  MAU  data  acquisiton-and- 
processing  system  are  shown  in  Fig.  3-1.  The  digitized  LDV  intensity 
versus  frequency  signal  along  with  its  coordinates  in  space  is  fed  into  the 
SEL  810  minicomputer.  Preprocessing  of  the  LDV  signal  is  carried  out  on 
the  minicomputer  utilizing  on-line  computer  programs  written  in  SEL  ma- 
chine language.  Information  from  the  SEL.  8 10  is  stored  on  magnetic  tape 
and  is  used  as  an  input  to  the  off-line  processing  algorithms.  Off-line  proc- 
essing of  the  LDV  signal  is  carried  out  on  a Univac  1108  computer  with 
programs  written  in  FORTRAN  language  and  using  card  inputs  with  informa- 
tion from  the  data  logs  to  supplement  the  data.  The  final  output  consists  of 
printouts  and  plots.  A description  of  the  data  logger  program  and  the  VAD 
program  and  the  operational  characteristics  of  these  programs  is  given  in 
the  following  sections. 

3.1  DESCRIPTION  OF  LASER  DOPPLER  VELOCIMETER  SOFTWARE  SYSTEM 

Data  acquisition  in  the  MAU  is  carried  out  by  the  SEL  Data  Logger 

program.  The  Data  Logger  program  preprocesses  and  records  the  LDV 

signal.  A flowchart  of  the  Data  Logger  program  is  given  in  Fig.  3-2.  A 

sweeping  spectrum  analyzer  is  used  to  detect  the  Doppler  shift  frequency. 

A diagram  of  the  output  of  the  spectrum  analyzer  is  shown  in  Fig.  3-3.  For 

each  sweep  of  10-,  20-,  or  50-millisecond  duration,  the  Data  Logger  saves  the 

maximum  amplitude  LDV  signal,  Ims . and  if®  corresponding  frequency,  Vmg, 

which  are  above  both  the  amplitude  and  frequency  thresholds.  The  definition 

of  1 and  V and  the  shape  of  the  characteristic  LDV  spectrum  are  shown 
ms  ms 

in  Fig.  3-3.  The  velocity  at  maximum  signal  intensity  is  V^g  and  is  taken 
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FIGURE  3-1.  GENERAL  ELEMENTS  OF  MOBILE  ATMOSPHERIC  UNIT 
DATA  ACQUISITION -AND-PROCESSING  SYSTEM. 


V = Line-of  - s lght  velocity  corresponding  to  highest 

” Doppler  frequency  shift  detectable  above  ampli- 

tude threshold 

V Line -of- s ight  velocity  corresponding  to  Doppler 
frequency  shift  at  maximum  signal  intensity 


Dopplf  F requency  Shift,  Af  (MHz) 


FIGURE  3-3.  TYPICAL  SPECTRUM  ANALYZER  OUTPUT  IN  VELOCITY 
AZIMUTH  DISPLAY  SCAN. 
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to  be  the  characteristic  line-of-sight  velocity  associated  with  the  flow  phe- 
nomenon. This  choice  is  made  because  the  maximum  signal  intensity  is 
obtained  from  the  center  of  the  focal  volume. 

The  output  from  the  Data  Logger  program  consists  of  vmg  as  a func- 
tion of  time  and  space.  From  the  output  of  the  Data  Logger,  the  wind  field 
can  be  reconstructed  using  off-line  processing  routines. 

Final  processing  of  the  LDV  measurements  is  carried  out  by  the  VAD 
program.  A macro  flowchart  of  the  VAD  program  is  shown  in  Fig.  3-4.  In 
this  off-line  program  the  array  of  vmS  values  which  is  a function  of  time 
and  space  is  processed  to  yield  the  three-dimensional  wind  field  in  the  VAD 
mode.  The  program  is  formulated  to  calculate  winds  for  both  the  translated 
and  non-translated  LDV  signal.  A translated  signal  is  provided  when  the  LDV 
system  includes  a frequency  translator  which  distinguishes  between  positive 
and  negative  values  of  line -of-s  ight  velocity.  A non-translated  signal  provides 
only  the  absolute  value  of  line-of-sight  velocity  as  shown  in  Fig.  2-13.  How- 
ever, during  the  course  of  this  research  effort,  all  of  the  data  acquisition- 
and-processing  were  done  in  the  non-translate  mode. 

Three  techniques  have  been  implemented  to  compute  the  three- 
dimensional  wind  components:  (1)  a peak  algorithm  where  the  magnitude  and 
location  of  the  peak  signal  in  the  sinusoidal  LDV  VAD  signature  are  used  to 
compute  the  velocity  components;  (2)  a spectral  processing  for  the  winds 
using  the  derectified  signal;  and  (3)  a sine  curve  fit.  The  final  output  is  a 
printout  (and  selected  plots)  of  the  wind  velocity  components  as  a function 
of  altitude  and  time.  Wind  velocity  components  are  given  for  both  rectilinear 
orthogonal  components  and  cylindrical  components. 

3.2  OPERATION  OF  LASER  DOPPLER  VELOCIMETER  SOFTWARE  SYSTEM 

Operation  of  the  Lockheed-Huntsville  MAU  involves  initialization  of  the 
SEL  Data  Logger  program  and  the  recording  of  the  LDV  signatures.  After 
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FIGURE  3-4.  MACRO  FLOWCHART  OF  OFF-LINE  VELOCITY  AZIMUTH 
DISPLAY  PROGRAM. 
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the  measurements  have  been  recorded  at  the  proper  threshold  settings,  the 
VAD  program  is  used  to  process  the  data  into  VAD  plots.  It  is  useful  to  ex- 
amine the  data-proces sing  operations  involved  in  determining  the  three- 
dimensional  wind  velocity  field  from  the  LDV  measurements  in  terms  of 
on-line  preprocessing  and  off-line  processing. 

3.2.1  On-Line  Data  Processing 

The  conical  scan  VAD  measurements  are  preprocessed  on  the  SEL 
computer  and  recorded  on  magnetic  tape.  The  output  from  the  SEL  computer 
consists  of  the  basic  Vmg  signal  as  well  as  additional  test  parameters  which 
are  listed  in  Table  3-1.  It  is  useful  to  consider  the  type  of  data  recorded  by 
the  SEL  during  normal  operations  in  the  VAD  mode. 


A dump  of  a sample  output  tape  from  the  Data  Logger  operating  in  the 
VAD  mode  is  shown  in  Fig.  3-5,  each  row  of  data  corresponding  to  information 
recorded  for  each  spectrum -analyzer  sweep.  The  information,  packed  into 
four  computer  words  on  tape,  is  separated  into  13  columns  in  the  printout  with 
interpretation  of  the  various  columns  being  as  follows: 


Column 


Interpretation 


1 (PEAK/MAX)  An  indication  as  to  whether  the  peak  Doppler  frequency 

(fpk  corresponding  to  Vpk  in  Fig.  3-4)  or  maximum 
signal  Doppler  frequency  (fms  corresponding  to  Vms  in 
Fig.  3-4)  is  used.  Indication  is  0 if  maximum  signal 
frequency  is  used  and  is  1 if  peak  velocity  signal  is 
used. 

2 (0  = 0)  A conical-scan  azimuth  reference  which  is  nonzero  when 

the  reference  switch  is  activated. 


3  (NOR/VAD)  An  indication  as  to  whether  the  system  is  operating  in 

the  normal  (vortex  track)  (1)  or  VAD  (0)  mode. 


4  (SWEEP  SPEED)  Sweep  speed  of  the  spectrum  analyzer  trace  in  msec/cm 


5 (DIG  TRCK) 


The  computer-calculated  estimate  of  f , (percent  of 

pK 


full  scale). 


6 (DATA  ACCEPT)  An  indication  that  there  is  (1)  or  is  not  (0)  an  output 

above  the  frequency  and  amplitude  thresholds  during  a 
sweep. 
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TABLE  3-1 

INPUT  FROM  SYSTEMS  ENGINEERING  LABORATORIES  COMPUTER 

Each  VAD  run  is  in  a separate  file.  The  following  data  are  recorded 
for  each  spectrum  analyzer  sweep. 

1.  Spectral  Sample  Count  Across  Spectrum  Analyzer  Sweep, 
Corresponding  to 

2.  Amplitude  at  Above  Point 

3.  Data  Acceptable  Flag 

4.  Flag  for  Spectrum  Analyzer  Sweep  Speed 

5.  Flag  for  Translator 

6.  Flag  for  Positive  or  Negative  Frequency  (used  only  when 
translator  is  used) 

7.  Flag  for  V . or  V 

a pk  ms 

8.  Flag  for  Conical  Scan  or  Normal  Scan 
For  Conical  Scan 

9.  Height  above  Van 

10.  Flag  for  Azimuth  Switch 

1 1 . Cone  Angle 


FIGURE  3-5.  DUMP  OF  SAMPLE  OUTPUT  TAPE  FROM  DATA  LOGGER  (With  System  Operating  i 
Velocity  Azimuth  Display  Mode). 


Column 


Interpretation 


7 (+/-) 


8 (TRANS/NON- 
TRANS) 


9 (SPECTRUM 
INTENSITY) 


An  indication  as  to  the  sense  of  the  Doppler  shift: 

1 r moving  toward  MAU,  0 - moving  away  from  MAU. 

An  indication  as  to  whether  or  not  a frequency  transla- 
tor was  incorporated.  During  this  investigation,  it  was 
not  incorporated.  1 = Yes,  0 = No. 

Peak  amplitude  of  the  Doppler  spectrum  in  region  above 
a frequency  threshold  (arbitrary  units  ranging  from  0 
to  1024). 


10  (NBR  ROTN) 

> l (ALT  (M)) 

12  (LTRNC  TRCK) 


Number  of  successive  VAD  scans  for  a particular  altitude. 

Altitude  of  VAD  for  particular  sweep  in  meters. 

On-line  frequency-tracker  estimate  of  f , (should  be 
approximately  equal  to  column  5).  P 


13  (CONE  ANGLE)  Half-angle  of  VAD  cone  in  degrees. 


3.2.2  Off-Line  Data  Processing 

Three  algorithms  for  calculating  the  mean  wind  speed  and  direction  from 
the  VAD  signature  have  been  developed.  For  each  of  these  algorithms,  mean 
wind  speed  and  direction  are  calculated  for  each  5-second  VAD  sweep.  Stan- 
dard deviations  of  wind  speed  and  direction  can  be  calculated  from  multiple 
VAD  sweeps.  The  data  output  of  the  LDV  system  operating  in  the  VAD  mode 
are  line-of-sight  velocities  measured  at  a selected  number  (usually  350  in 
the  current  Lockheed  system)  of  distinct  points  around  the  VAD  cone.  Recall 
that  the  line-of-sight  velocity  signature  is  theoretically  sinusoidal  in  the  VAD 
mode  (cf.,  Fig.  2-12). 


For  all  the  processing  algorithms,  preprocessing  of  the  data  occurs; 
the  preprocessing  includes: 

a.  Save  line-of-sight  velocities  for  one  rotation  of  scanner. 

b.  If  two  or  more  rotations  occur  at  the  same  altitude,  average 
with  previous  rotation(s). 

c . Assign  azimut-i  angle  to  each  point  (assuming  constant  rotation 
rate). 
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d.  Edit  data  to  eliminate  spurious  points.  Eliminated  points  are  de- 
clared unacceptable. 

e.  Plot  line-of-s ight  velocity  versus  azimuth  angle. 

The  edit  criterion  for  the  elimination  of  spurious  points  is  that  the  i**1 
point  is  eliminated  if 


and 


v . - v . . , 
r . i r , t + 1 


> .2  v 


r , i + 1 * 


v . - v . 
r , l r , i - 1 


> .2  v 


r, i-  1 ’ 


A sample  plot  of  unedited  lin e -of - s ight  velocity  versus  azimuth  angle  is  shown 
in  Fig.  3-6.  Unacceptable  data  points  are  shown  as  zero  velocity. 


3. 2. 2.1  Peak  Algorithm 

For  the  calculation  of  wind  velocity  by  the  peak  algorithm,  the  pro- 
cedure is: 


, that  occur 


a.  Filter  data  with  an  n-point  moving  average. 

b.  Identify  the  two  peak  velocity  points,  V j and  V 

it  a minimum  of  90  degrees  apart.  ^ E 

c.  Compute  the  magnitude  of  the  horizontal  component  of  the  wind 

V , + V , 

V E] sA  . 

h 2 sina 


d.  Compute  horizontal  wind  angle  with  help  of  estimated  wind 
direction. 

e.  Compute  the  vertical  component  of  the  wind 

V . - V , 

w = _Ei El  . 

v 2 cosa 

f.  Dorectify  VAD  signal  if  no  translator  is  present  and  plot  derectified 
signal.  The  signal  is  derectified,  so  that  the  positive  peak  of  the  derectified 
signature  Is  the  peak  which  is  closer  to  the  estimated  wind  direction  recorded 
by  the  operator  at  the  time  the  data  are  measured. 


MRitacr*.  MU  m 


FIGURE  3-6.  SAMPLE  OUTPUT  PLOT  FROM  VELOCITY  AZIMUTH  DISPLAY 
PROGRAM  (Unedited). 


The  purpose  of  the  filter  is  the  elimination  of  nonuniformities  in  the  data. 
This  is  especially  important  for  the  peak  algorithm.  The  individual  data  points 
are  filtered  with  an  n-point  moving  average  filter  (n  is  usually  21).  Thus,  each 
line -of-sight  point  is  filtered  by 


r,  i 


n-  1 


+ V . + V .+V  + 

r,  t-1  r,  l r , i + I 


+ V 


r,  i + 


n - 1 

2 J 


/n. 


where  V . is  the  filtered  value  of  the  line-of-sight  velocity  to  be  used  in  further 
r 1 1 

calculations.  A plot  of  the  filtered  line-of-sight  velocities  for  a 21-point  filter 
is  shown  in  Fig.  3-7.  Additional  samples  of  the  LDV  signature  (including  raw 
data,  filtered  data,  and  derectified  data)  are  presented  in  Appendix  A. 

When  the  LDV  data  are  measured,  an  approximation  of  the  wind  angle  is 
recorded.  The  calculated  wind  angle  is  the  azimuth  angle  of  the  peak  which  is 
Closer  to  the  estimated  wind  angle.  For  small  values  of  vertical  velocity  the 
wind  angle  plus  90  deg  is  the  angle  at  which  the  line-of-sight  velocity  is  theo- 
retically zero.  This  angle  is  used  for  the  derectification  of  the  line-of-sight 
signal.  A plot  of  the  derectified  (edited,  but  unfilte.-ed)  line-of-sight  velocity 
is  presented  in  Fig.  3-8. 

3. 2. 2. 2 Fourier  Coefficient  Algorithm 

The  Fourier  coefficient  algorithm  (or  spectral  algorithm)  computes  the 
fundamental  harmonic  of  the  line-of-sight  velocity.  The  Fourier  series  for  a 
generalized  periodic  function  is 


f(x) 


ao 

Z(a  cos  nx 
n 

n=  1 


+ b^  sin  nx  ) , 
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FIGURE  3-7.  FILTERED  LINE-OF-SIGHT  VELOCITY  FOR  VELOCITY 
AZIMUTH  DISPLAY  MODE. 


>« 


where 


a„  = */ 


/in 

v' 


(0)  cos  n6  d0  , 


bn  = 1/x 


rin 

J v< 


(0)  sin  nfl  d6. 


The  integrals  are  performed  by  trapezoidal  integration  on  the  discrete 
points  of  the  edited  and  derectified  velocity  signature  (cf.  Fig.  3-8).  The 
calculation  procedure  for  the  spectral  algorithm  is: 

a.  Compute  tne  Fourier  coefficients  for  the  dc  level  (i.e., 
n=0)  and  the  first  three  harmonics  (n  = 1,  2,  3)  using  the 
above  equations. 

b.  Compute  the  magnitude  of  the  component  of  wind  velocity 

\A\ + 

AT  — A 1 


c.  Compute  the  horizontal  wind  angle 


angle^  = Atan(bj/aj). 


d.  Compute  the  vertical  component  of  wind  velocity 


2 cosa 


3. 2.2. 3 Sine-Curve-Fit  Technique 


The  line-of-sight  velocity  signature  is  sinusoidal  for  a uniform  wind. 
Therefore,  the  best  sinusoidal  wave  which  fits  the  data  in  a least-squares 
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sense  is  determined  by  finding  the  values  of  the  coefficients  A,  B, 
and  C which  minimize  (summing  over  all  acceptable  data  points) 


^ * (V.  - C - A cosf^  - B sinS^)"  , 
i 

where  V is  the  line -of  - sight  velocity  (derectifii  ) at  point  i,  and  8^  is  the 
azimuth  at  point  i.  A,B,  and  C are  obtained  from 


f^Pcos^O^A  + f y *|cos8.  sinB^B  + ^y^cosB  jc 
j^Scos8.  sinO^A  + [^sin^B  + l?  sinft^C 
[X>ei]A  -t^  ^sinS^B  4-  nC 


Y!  V.  cos8if 


EV.  sin8., 
I I 

2X- 


The  steps  for  calculating  the  wind  using  the  least-squares  algorithm  are 


a.  Find  the  least-squares  curve  fit  for  a sine  wave  to  minimize 


£(Vt  - C - A cos6^  - B sin6^  )^, 


wher  e 

is  line-of- sight  velocity  at  point  i, 
is  azimuth  at  point  i,  and 
C,  A,  and  B are  coefficients  to  be  calculated. 


b.  Compute  the  magnitude  of  the  horizontal  component  of  velocity: 


V 


h 


sin  a 


c.  Compute  the  horizontal  angle: 

Angle^  = Atan(B/A). 
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d.  Compute  the  vertical  component  of  wind  velocity: 


w 

v 


-C 

c os  a 
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4.  DATA  COLLECTION  AT  NATIONAL  SEVERE  STORMS 
LABORATORY  TEST  SITE 

4.1  TEST  DESCRIPTION 

During  June  1976,  the  mobile  laser  Doppler  velocirr r.ter  was  deployed 
adjacent  to  the  WKY-TV  tower  north  of  Oklahoma  City.  The  meteorological 
instruments  on  the  481-m  tower  are  operated  by  the  National  Severe  Storms 
Laboratory  (NSSL)  of  the  National  Oceanic  and  Atmospheric  Administration. 

The  tower  was  instrumented  with  propeller  anemometers  at  the  surface  and 
at  altitudes  of  25,  45,  89,  177,  266,  and  444  m. 

A diagram  of  the  test  site  is  shown  in  Fig.  4-1,  and  a table  of  test  runs 
is  shown  in  Table  4-1.  Meteorological  data  from  the  tower  were  recorded  at 
2-sec  intervals.  Table  4-2  summarize  the  data  recorded.  The  data  recorded 
at  2-sec  intervals  were  averaged  over  appropriate  time  intervals  for  com- 
parison with  LDV-measured  wind. 

By  standard  meteorological  convention,  the  wind  direction  is  the  di- 
rection from  which  the  wind  comes.  Therefore,  when  wind  is  given  in  speed- 
direction  coordinates,  the  direction  is  the  direction  from  wh  'i  the  wind  comes, 
measured  from  true  north.  The  coordinates  chosen  for  the  test  are  shown  in 
Fig.  4-1.  For  wind  expressed  in  rectangular  coordinates  the  positive  u com- 
ponent represents  wind  coming  from  the  positive  x direction  (i.e.,  from  north 
to  south).  Similarly,  a positive  v component  represents  wind  coming  from  the 
positive  y component  (i.e.,  from  west  to  oast).  For  example,  a wind  coming 
from  045  deg  has  a positive  u component  and  a negative  v component.  The 
vertical  component,  w,  is  positive  for  upward  air  motion. 

4.2  TEST  OBJECTIVES  AND  OPERATING  PROCEDURE 

The  primary  objective  of  the  test  was  the  verification  of  the  capability 
of  the  laser  Doppler  velocimeter  to  perform  accurate  wind  measurements  to 
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ocation  of  Lockheed 


TABLE  4-1.  SUMMARY  OF  LASER  DOPPLER  VELOC1METER  WIND  MEASURE- 
MENTS AT  NATIONAL  SEVERE  STORMS  LABORATORY  TEST  FACILITY. 


Date 

Run  No. 

Altitude 

Start  Time 

Stop  Time 

Comments 

AGL 

(CST) 

(CST) 



WJMJMim 

266 

16:22 

BTFTH 

6/8/76 

100 

444 

09:25 

09:54 

201 

444 

10:01 

302 

266 

1 1:22 

303 

177 

12:54 

304 

45 

14:37 

16:21 

100 

45 

11:14 

12:10 

201 

45 

LU27 

12:45 

302 

177  _J 

12:58  _ 

13:37 

100 

444 

22:22 

24:27 

Printed  Time  Wrong 

2 Altitudes 

201 

266 

24:27 

02:27 

302 

177 

02:30 

03:32 

Lost  (Bad  Tape) 

403 

45 

03:38 

04:37 

Lost  (Bad  TaDe) 

6/1 1/76 

100 

45 

17:55 

hi?  in 

201 

177 

19:05 

■2  is 

302 

266 

20:40 

403 

444 

22:47 

00:25 

504 

444 

00:31 

01:17 

Printed  Time  Wrong 

2 Altitudes 

6/13/76 

100 

444 

17:53 

18:49 

201 

444 

18:55 

18:59 

6/ 14/76 

100 

; H 

11:15 

■fm 

201 

12:37 

302 

-if'l 

■mTfH 

6/14/76 

100 

444 

21:45 

22:57 

Lost  (Bad  Tape) 

201 

444 

23:10 

23:23 

Lost  (Bad  Tape) 

WJMEMBa 

100 

45 

01:55 

DXSXj 

100 

45 

09:35 

■B52SI 

V 

201 

177 

11:19 

302 

266 

13:03 

led 

' 

403 

444 

14:47 

15:26 

J|:1 

101 

444 

19:26 

22:41 

Wind  Shifts  at  21:42 
and  22: 15 

202 

444 

22:44 

23:07 

303 

444 

23:09 

23:15 

23:18 

6/21/76 

101 

. 

444 

08:39 

09:44 

202 

444 

09:46 

10:32 

303 

266 

10:36 

12:07 

404 

177 

12:16 

13:38 

505 

45 

VTttTH 

16:14 

6/23/76 

101 

444 

21:13 

21:53 

Full  Spectrum  Data 

202 

444 

21:53 

22:06 

Rain  Began 

303 

444 

22:38 

22:45 

404 

266 

22:46 

23:41 

505 

45 

00:14 

00:53 

Rain  Stopped 
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TABLE  4-2.  METEOROLOGICAL  DATA  RECORDED 
FROM  WKY-TV  TOWER 


Parameter 

Unit 

Altitude  (m) 

Wind  Speed 

, . * 

Wind  Direction 

Temperature 

Wet  Bulb  Temperature 

Relative  Humidity 

** 

Vertical  Velocity  Component 

Pressure 

m/  sec 

deg 

deg  C 

deg  C 

percent 

m/  sec 

millibar 

Surface,  26,  45,  89,  177,  266,  444 

Surface,  26,  45,  89,  177,  266,  444 

Surface,  26,  45,  89,  177,  2 66,  444 

Surface 

89, 266,  444 

26,  45,  89,  177,  266,  444 

Surface,  444 

*Wind  speed  and  direction  were  measured  by  Bendix  Model  120  aerovanes. 

**Vertical  component  of  velocity  was  measured  by  R.M.  Young  Model  1200 
propeller  anemometers. 
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an  altitude  of  444  m.  As  shown  in  Table  4-1,  measurements  were  made  at 
45,  177,  <166,  and  444  m.  For  all  runs,  measurements  were  made  at  one  alti- 
tude at  a time,  although  the  altitude  was  occasionally  changed  during  a run. 
Therefore,  all  averaged  data  are  based  on  12  scans  per  minute.  A secondary 
objective  was  the  measurement  of  wind  phenomena  in  thunderstorms  and  a 
demonstration  of  the  ability  of  the  LDV  to  measure  wind  in  rain.  A small 
quantity  of  data  was  obtained  in  rain. 

4.3  ANALYSIS  OF  WIND  MEASUREMENTS 
4-3.1  Effect  of  Measurement  Altitude 


Appendix  A shows  several  typical  VAD  signatures  similar  to  those  shown 
in  Figs.  3-6,  3-7,  and  3-8.  Figures  4-2  and  4-3  show  signatures  (unedited  and 
unfiltered)  for  altitudes  of  45  and  444  m.  Intuitively,  it  might  be  expected  that 
the  data  from  444  m would  exhibit  more  scatter  than  the  data  measured  at  45  m 
It  is  noted  that  the  focal  volume  (i.e.,  the  volume  of  space  from  which  meas- 
urements are  taken)  includes  altitudes  of  45  m +0.8  m for  the  45-m  data,  and 
includes  altitudes  of  444  m _+82  m for  the  444-m  data.  However,  the  signatures 
show  that  the  sinusoidal  curve  for  the  high  altitude  is  very  well  defined,  and  no 
ambiguity  about  the  data  occurs.  It  is  noted  that  wind  variation  with  altitude  is 
much  less  significant  at  the  higher  altitude  (i.e.,  out  of  the  earth's  boundary 
layer)  and  this  may  account  for  the  well-defined  curve  of  Fig.  4-3. 

4.3.2  Presentation  of  Measurements  Made 

The  winds  measured  by  the  LDV  are  tabulated  for  typical  runs  in 
Appendix  B.  A sample  of  the  output  is  shown  in  Table  4-3. 

The  date  is  the  test  site  date  on  which  the  run  was  initiated.  All  times 
are  Central  Standard  Time.  The  height  is  the  indicated  altitude  (m)  on  the 
laser  focusing  system.  Before  each  test  run,  the  laser  ranging  system  is 
calibrated  to  assure  that  the  actual  altitude  of  the  data  is  the  desired  altitude. 
The  following  interpretations  are  placed  on  the  columns  in  Table  4-2. 
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FIGURE  4- 


. LINE-OF-SIGHT  VELOCITY  SIGNATURE  AT  45-M  ALTITUDE 
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FIGURE  4-3.  LINE-OF -SIGHT  VELOCITY  SIGNATURE  AT  444-M  ALTITUDE. 
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TABLE  4-3.  SAMPLE  TABULATED  WIND  DATA 
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TABLE  4-3  (Concluded) 
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4.3.2. 1 One-Minute  Means 


One-minute  means  were  calculated  for  each  minute  of  the  15-minute 
time  period.  The  minutes  of  the  time  period  are  tabulated  in  the  first 
column.  Each  of  the  numbers  in  this  group  is  an  average  of  the  12  VAD 
scans  that  occur  at  the  given  altitude  during  1 min. 

4. 3.2. 1.1  Peaks 

The  1-min  mean  u,  v,  and  w components  of  wind  as  calculated  by  the 
peak  algorithm  are  tabulated.  The  w component  is  the  vertical  component 
and  is  positive  vertically  upward.  Speeds  are  given  in  centimeters  per 
second  and  angles  are  given  in  degrees. 

The  data  calculated  from  the  basic  VAD  data  are  horizontal  wind  speed, 
wind  direction,  and  vertical  component  of  wind  velocity.  The  u and  v com- 
ponents are  calculated  from  horizontal  wind  speed  and  direction.  The  speed 
is  the  1-min  average  of  the  horizontal  component  of  velocity,  and  "TH"  is  0q, 
the  direction  of  the  horizontal  component  of  wind.  The  u and  v averages  are 
calculated  by  averaging  the  u and  v values  of  the  individual  VAD  scans  rather 
than  from  averaged  values  of  u and  v.  Thus,  if  S is  wind  speed, 

ui  = - si  sin9o,i  • 

and 

v.  = - S.  cosfl.  . , 

l l 0,  i 

where  the  i subscript  refers  to  the  individual  VAD  scans.  The  averages  for 
n scans  in  the  1-min  average  are  (for  data  discussed  herein,  n=  12) 

S = £ S./n  , n) 

»o  = Ev/“' 

1 = 2 V"  • 

and 

v = £ v./n  . 
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It  is  noted  in  general  that 


u t - S sin0Q  , 
v / - S cos0Q  , 

and  

S / Vu2  + v2  . 

The  averaging  characteristic  represented  by  the  above  equations  also 
results  from  resolving  cup  and  vane  anemometer  data  into  components  and 
averaging  or  when  comparing  averaged  data  from  cup  and  vane  anemometers 
with  averaged  data  from  propeller  anemometers.  The  characteristic  is  par- 
ticularly pronounced  in  light  and  variable  wind  conditions. 

4.3.2. 1.2  Fourier  Coefficients 

Results  are  shown  for  wind  calculations  performed  using  the  spectral 
algorithm.  The  symbols  have  the  same  meaning  as  for  the  peak  algorithm. 
The  " 2D"  is  the  magnitude  of  the  second  harmonic  and  is  the  average  of 
Vaf  + b^  where  a,  is  the  second  Fourier  cosine  coefficient  and  b,  is  the 

“*  * * i > 2 

second  Fourier  sine  coefficient.  Similarly,  "3D"  is  the  average  of  Va^  +•  b^. 
The  magnitude  of  these  numbers  is  an  indication  of  turbulence  with  the  fre- 
quency of  the  turbulence  increasing  with  the  number  of  the  harmonic. 

4. 3.2. 1.3  Sine  Fit 

Results  are  shown  for  the  least -squares  sine  fit  algorithm  with  similar 
meanings  of  the  symbols  as  for  the  other  two  algorithms.  The  "SP"  is  the 
standard  deviation  of  the  actual  line -of-sight  velocity  points  referenced  to  the 
best  sine  fit.  As  such,  it  is  an  indication  of  atmospheric  turbulence. 

4.  3. 2. 2 Cumulative  Means 

The  cumulative  means  are  averages  of  individual  VAD  scans  taken 
over  the  time  period  indicated.  For  example,  the  3-min  cumulative  mean 
is  the  average  of  the  VAD  scans  taken  over  the  first  3 min  of  the  15-min 
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time  period.  They  are  not  averages  of  the  1-min  means.  The  two  averaging 
methods  would  be  identical  if  the  number  of  VAD  scans  in  each  1-min  average 
were  always  identical.  The  15-min  cumulative  means  are  the  averages  over 
the  15-min  averaging  period. 

4. 3. 2. 3 1 -Min  Standard  Deviations 

This  section  lists  the  standard  deviation  of  each  of  the  variables  for 
each  1-min  time  period.  Each  standard  deviation  is  the  standard  deviation  of 
the  six  VAD  scans  in  the  1-min  period. 

4. 3. 2 .4  Cumulative  Standard  Deviations 

The  cumulative  standard  deviations  are  analogous  to  the  cumulative 
means.  They  are  the  standard  deviations  of  the  individual  VAD  scans  taken 
over  the  time  period  indicated.  In  general,  the  cumulative  averages  of  the 
standard  deviations  are  not  averages  of  the  1-min  means  because  of  the 
mathematical  definition  of  standard  deviation. 

4.3.3  Averaging  Time  Selection 

In  the  tabulated  data,  the  15-min  averaging  periods  were  chosen  to  start 
on  the  quarter  hours.  However,  if  the  start  and  stop  time  of  an  individual  run 
were  such  that  more  15-min  averages  could  be  extracted  from  the  available 
data  by  starting  the  averaging  periods  at  times  other  than  the  quarter  hour, 
the  starting  time  which  yielded  the  most  15-min  averages  was  used.  Data 
runs  lasting  less  than  15  min  were  not  processed. 

4.4  DATA  ANALYSIS 

4.4.1  Comparison  of  Results  for  Various  Algorithms 

From  the  tabulated  data,  it  is  noted  that  the  spectral  algorithm  gives  a 
value  of  wind  speed  which  is  higher  than  that  for  the  other  two  algorithms. 

In  previous  studies  (e.g..  Ref . 6),  the  spectral  algorithm  gave  a lower  value 
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« 


of  wind  speed.  For  the  previous  studies,  the  Fourier  coefficients  were  c 'cu- 

» 

lated  by  summation  around  the  VAD  scan  with  an  application  of  an  analytically- 
derived  correction  factor  to  compensate  for  the  absence  of  data  points  which 
lie  below  the  velocity  threshold.  In  general,  the  wind  speeds  obtained  by  this 
method  were  lower  than  the  wind  speeds  measured  by  the  other  two  algorithms 
as  well  as  those  measured  by  the  tower. 

In  order  to  correct  this  condition,  the  method  of  calculating  the  Fourier 
coefficients  was  changed  to  trapezoidal  integration  of  the  edited,  derectified 
signal.  This  technique  obviates  the  need  for  a correction  factor  because  it 
effectively  places  a straight  line  across  the  gap  in  the  derectified  signature 
caused  by  the  velocity  threshold  (cf..  Fig.  3-8).  However,  the  technique  causes 
the  wind  speed  to  be  biased  high.  The  cause  of  the  bias  is  yet  to  be  identified. 

A comparison  between  laser-measured  winds  and  tower-measured  winds 
is  shown  for  the  three  algorithms  in  Figs.  4-4,  4-5,  and  4-6.  The  peak  algo- 
rithm exhibits  a slightly  high  bias,  the  spectral  algorithm  exhibits  a significantly 
high  bias,  and  the  least-squares  sine  algorithm  exhibits  little  systematic  bias. 

It  is  noted  that  two  significant  wind  shifts  occurred  during  the  run,  and  a wide 
range  of  wind  speeds  is  shown.  The  comparison  presented  in  Ref.  6 showed 
that  the  sine  algorithm  appears  to  be  the  most  reliable  data  processing  algo- 
rithm. This  conclusion  is  confirmed  by  Figs.  4-4,  4-5,  and  4-6.  The  con- 
clusion of  a high  bias  for  the  modified  Fourier  algorithm  (i.e.,  trapezoidal 
integration)  is  based  on  comparison  with  the  sine  algorithm,  which  had  pre- 
viously been  validated  in  Ref.  6,  as  well  as  comparison  with  anemometer- 
measured  data.  Since  the  sine  algorithm  appears  to  be  the  most  reliable  data 
processing  algorithm,  further  comparison  with  the  meteorological  tower  is  pre- 
sented for  the  sine  algorithm  only. 

4.4.2  Effect  of  Averaging  Period  on  Wind  Results 

■ The  convergence  of  the  wind  data  with  increasing  averaging  times  is  shown 

in  Figs.  4-6  through  4-9.  Similar  data  for  selected  runs  for  1-,  3-,  6-,  9-,  12-, 

. and  15-min  averages  are  shown  in  Appendix  C.  Wind  direction  data  are  also 

presented.  The  scatter  for  wind  direction  is  significantly  less  than  that  for 
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FIGURE  4-4.  1 - MIN  AVERAGES  OF  WIND  SPEED  USING  PEAK  ALGORITHM. 
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FIGURE  4-5.  1 -MIN  AVERAGES  OF  WIND  SPEED  USING  SPECTRAL 

ALGORITHM. 
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FIGURE  4-6.  1 - MIN  AVERAGES  OF  WIND  SPEED  USING  SINE 

ALGORITHM. 
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FIGURE  4~7.  3-MIN  AVERAGES  OF  WIND  SPEED  USING  SINE  ALGORITHM. 
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wind  speed.  Appendix  D shows  3-min  and  15-min  average  comparisons  for  runs 
not  presented  in  Appendix  C. 

4.4.3  Comparison  Between  Laser -Measured  Winds  and  Anemometer- 
Measured  Winds 

A review  of  the  data  shown  in  Appendixes  C and  D shows  a significant  differ- 
ence between  the  wind  speed  measured  by  the  LDV  and  that  measured  by  the  tower- 
mounted  anemometers.  Wind  directions  show  good  agreement.  Table  4-4  shows 
a synopsis  of  the  data  in  Appendixes  C and  D,  including  the  deviation  of  the  laser- 
measured  wind.  Several  characteristics  of  the  data  are  noted: 

The  deviation  does  not  seem  to  be  dependent  on  wind  speed.  Thus,  the 
deviation  appears  to  be  a bias,  rather  than  a proportionality  problem. 

The  deviation  appears  to  be  more  significant  for  lower  altitudes  than 
for  the  444-m  altitude. 

Most  of  the  data  are  for  wind  directions  of  150  to  200  deg.  The  data 
for  wind  azimuth  greater  than  200  deg  show  no  systematic  deviation, 
suggesting  that  the  deviation  may  be  directionally  dependent. 

The  source  of  the  deviation  cannot  be  positively  identified.  However, 
reasons  for  the  deviation  can  be  postulated  as  (in  decreasing  order  of  likelihood): 
inherent  errors  in  the  wind  measurements  made  by  the  WKY-TV  tower  ane- 
mometers, anemometer  errors,  local  terrain  effects,  altitude  calibration  errors 
of  the  LDV,  LDV  cone  angle  errors,  and  errors  in  the  LDV.  Some  of  these  pos- 
sible sources  for  the  discrepancy  have  previously  been  reported  for  the  WKY-TV 
tower. 

4.4.3. 1 Inherent  Errors  in  Tower-Measured  Winds 

The  most  likely  source  of  the  discrepancy  is  a tower  characteristic  which 
causes  the  wind  measured  at  the  tower  to  be  slightly  less  than  that  measured  in 
the  vicinity  of  the  tower.  It  is  noted  that  previous  tests  of  the  LDV  at  other 
locations  have  failed  to  indicate  any  systematic  discrepancy  between  laser- 
measured  winds  and  the  anemometer-measured  winds,  whereas  a discrepancy 
between  anemometer-measured  winds  and  winds  measured  by  different  tech- 
niques has  previously  been  reported  for  the  WKY-TV  tower. 
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TABLE  4-4.  DEVIATION  BETWEEN  LASER -MEASURED  WIND  SPEED 
AND  ANEMOMETER-MEASURED  WIND  SPEED. 


A comparison  between  laser-measured  wind  speed  and  anemometer- 
measured  wind  speed  as  measured  at  Table  Mountain,  Colorado,  is  shown  in 
Fig.  4-10  (from  Ref.  6).  From  these  data,  it  is  seen  that  agreement  is  usually 
within  _f30  cm/sec  and  is  always  within  1.5  m/sec.  The  altitudes  roughtly  corre- 
spond with  the  lower  two  altitudes  of  the  WKY-TV  tower.  The  validity  of  the 
LDV  and  the  least-squares  sine  fit  algorithm  has  been  established  from  the 
Table  Mountain  test. 

Reference  9 reports  a test  in  which  the  winds  measured  by  anemometers 
on  the  WKY-TV  tower  were  compared  with  observations  from  double-theodolite 
pilot  ballons  and  tetroons.  The  results  from  Ref.  9 are  shown  in  Fig.  4-11. 

Several  observations  about  the  data  are  significant: 

The  wind  direction  is  approximately  equal  to  that  during  the  LDV  tests. 

The  discrepancy  between  the  anemometer-measured  wind  speed  and 
balloon-observed  wind  speed  ranges  from  approximate  agreement  to 
the  anemometers  underestimating  wind  by  approximately  2.5  m/sec. 

The  same  discrepancy  range  was  observed  with  the  LDV. 

These  results  suggest  that  the  discrepancy  between  laser-measured  winds  and 
anemometer-measured  winds  is  due  to  the  underestimating  of  winds  in  the 
vicinity  of  the  tower  by  the  tower-mounted  anemometers. 

4.4.  3.2  Anemometer  Errors 

Table  4-3  shows  that  the  greatest  deviation  between  the  laser-measured 
wind  and  the  anemometer-measured  wind  occurs  at  an  altitude  of  266  m.  The 
anemometer  data  also  show  an  anomaly  at  266  m.  Table  4-5  shows  the  1-min 
mean  wind  speeds  as  measured  by  the  anemometers  during  the  early  part  of 
run  13-3.  At  10:40,  the  measured  wind  at  266  m is  significantly  less  than  that 
at  177  m or  444  m.  Although  this  condition  does  not  occur  universally,  it  occurs 
for  a significant  number  of  1-min  averages  throughout  the  run.  The  magnitude 
of  the  deviation  and  its  frequency  of  occurrence  suggests  that  at  least  one  anomaly 
exists  in  the  anemometer  data.  The  sign  of  the  anomaly  at  166  m in  Table  4-5 
supports  the  fact  that  the  greatest  deviation  between  laser-measured  wind  and 
anemometer-measured  wind  occurs  at  that  altitude. 
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Laser  Wind  Speed  (m/sec) 


FIGURE  4-10.  COMPARISON  OF  LASER  AND  TOWER  15-MIN  MEAN  WIND 
SPEED  USING  LEAST-SQUARES  SINE  ALGORITHM. 
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FIGURE  4-11.  MEAN  VARIATION  OF  WIND  SPEED  (LEFT)  AND  DIRECTION 
(RIGHT)  WITH  HEIGHT  FOR  2-HR  INTERVALS,  AS  DETERMINED  FROM  AERO- 
VANES  3 M UPWIND  OF  THE  WKY  TOWER  (SOLID  LINES)  AND  DOUBLE- 
THEODOLITE PILOT-BALLOON  OBSERVATIONS  AT  6-MIN  INTERVALS 
WITHIN  2 KM  OF  THE  TOWER  (DASHED  LINES).  The  letters  G and  T indicate, 
respectively,  the  wind  determined  from  Gill  anemometers  3 m upwind  of  the  tower 
and  from  tetroons  passing  by  the  tower  (the  latter  an  average  for  the  40-min  tra- 
jectory segment  centered  on  the  tower  (from  Ref.  9). 
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4.4.3.  3 Local  Terrain  Effects 


Figure  4-1  shows  the  end  of  a gully  at  a bearing  of  approximately  160  deg 
from  the  tower.  Reference  8 shows  that  the  gully  extends  for  2.6  km  south- 
southeast  of  the  tower  and  has  in  elevation  of  approximately  30  m below  that 
of  the  tower  at  a distance  of  1.6  km  south-southeast  of  the  tower.  For  most  of 
the  test,  the  wind  was  blowing  up  the  gully.  The  terrain  effects  of  the  gully 
could  cause  the  wind  measured  at  a point  (i.e.,  the  tower)  to  be  different  from 
that  measured  over  an  area  (i.e.,  the  LDV).  This  possibility  is  supported  by 
the  fact  that  no  systematic  deviation  exists  when  the  wind  azimuth  is  greater 
than  200  deg.  However,  the  data  for  wind  azimuth  other  than  150  to  200  deg 
are  too  sparse  to  draw  meaningful  conclusions. 

The  possibility  of  local  terrain  effects  is  further  supported  by  the 
vertical  component  of  wind  as  measured  by  the  tower -mounted  anemometers. 
Typical  magnitudes  of  the  vertical  component  of  wind  at  the  26 - m altitude 
are  .7  to  .8  m/sec,  and  they  are  both  positive  and  negative.  Vertical  wind 
components  above  1 m/sec  are  not  uncommon,  particularly  at  the  higher 
altitudes.  The  magnitude  of  the  vertical  velocity  component  often  approxi- 
mates 10%  of  the  horizontal  wind  speed.  Table  B-l  shows  magnitudes  of 
the  vertical  velocity  component  which  are  significantly  greater  than  those 
measured  in  other  tests.  It  is  noted  that  when  the  wind  shifted  to  205  deg 
(page  B-  10),  the  vertical  component  became  insignificant.  This  suggests 
a local  terrain  effect. 

4.4. 3.4  Laser  Doppler  Velocimeter  Altitude  Calibration  Errors 

Before  every  run,  the  LDV  ranging  system  was  calibrated  to  assure 
that  the  actual  altitude  was  the  desired  altitude  of  measurement.  From  the 
data  shown  in  Table  4-5,  if  anomalous  points  are  removed,  the  variation  in 
wind  speed  with  change  of  altitude  is  not  sufficient  to  cause  deviations  as 
large  as  shown  in  Table  4-4,  even  if  altitude  calibration  was  wrong.  There- 
fore, LDV  altitude  calibration  errors  cannot  be  the  source  of  the  deviations 
shown  in  Table  4-4. 


66 


4.4.  3.5  Laser  Doppler  Velocimeter  Cone  Angle  Errors 


As  discussed  in  Section  3. 2. 2. 3,  the  equation  for  the  magnitude  of  the 
horizontal  component  of  wind  is 

v _ >IaZ  t-  BZ 

h ” sina  ’ 

where  a is  the  cone  angle  and  A and  B are  defined  in  Section  3.  Z.  2.  3.  Obviously, 
an  error  in  a will  cause  an  error  in  the  calculated  value  of  wind  speed  with  the 
error  being  proportional  to  the  wind  speed.  Since  the  cone  angle  is  constant 
from  day  to  day  (the  elevation  mirror  is  held  firm  by  a brake),  this  error  would 
not  show  variations  from  day  to  day.  However,  Table  4-4  shows  that  day-to- 
day  variations  in  the  wind-speed  deviation  occur.  Therefore,  it  is  concluded 
that  cone  angle  error  is  not  the  source  of  the  deviation  between  laser-measured 
wind  speed  and  anemometer-measured  wind  speed. 

4.4. 3.6  Laser  Doppler  Velocimeter  Errors 

It  is  possible  that  LDV  system  errors  and/or  errors  in  the  VAD  concept 
as  implemented  in  the  Lockheed  LDV  are  the  source  of  the  discrepancies  be- 
tween laser ‘measured  winds  and  anemometer-measured  winds.  However,  in 
three  other  laser-anemometer  wind  comparisons,  the  agreement  has  been  very 
good.  Reference  6 reports  a comparison  conducted  at  the  Table  Mountain  test 
site  of  the  National  Oceanic  and  Atmospheric  Administration  near  Boulder, 
Colorado.  Figure  4-10  shows  a sample  comparison  between  the  laser-measured 
wind  speed  and  the  anemomete r - mea sured  wind  speed.  The  discrepancy  is  much 
less  than  that  shown  in  Table  4-4  for  the  WKY-TV  comparisons.  Reference  10 
describes  a test  conducted  at  Otis  AFB,  Massachusetts,  in  September  1976. 

A comparison  between  laser-measured  wind  speed  and  anemometer-measured 
wind  speed  (Ref.  11)  showed  excellent  agreement.  A similar  test  in  1977  showed 
excellent  agreement  with  laser-measured  winds  calculated  in  real  time  (Ref.  11). 

4.4.3  Turbulence 


Figure  4-12  shows  the  VAD  signature  (21 -point  moving  average)  for  con- 
ditions which  were  extremely  turbulent.  Other  signatures  taken  for  this  run 
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FIGURE  4-12.  LINE-OF-SIGHT  VELOCITY  SIGNATURE  IN  HEAVY 
TURBULENCE. 
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show  the  same  characteristic.  The  fluctuations  may  be  identified  as  turbulence 
because  they  are  not  averaged  out  by  a 21-point  moving  average  (as  random  sys- 
tem noise  would  be),  the  other  VAD  signatures  measured  near  the  same  time 
exhibit  similar  characteristics,  and  the  signature  characteristics  do  not  occur 
as  significantly  for  other  runs. 

The  column  labeled  "SP"  on  the  tabulated  wind  data  (Table  4-3  and 
Appendix  B)  is  the  standard  deviation  of  the  difference  between  the  actual  line- 
of-sight  velocity  and  the  sine  wave  calculated  by  the  least-squares  curve  fit. 

As  such,  it  is  a measure  of  atmospheric  turbulence.  For  example,  for  the  1-min 
average  from  which  Fig.  4-12  was  taken,  SP  has  a value  of  73  cm/ sec.  For  less 
turbulent  data,  SP  has  a value  between  20  and  50  cm/sec.  Therefore,  the  value 
of  SP  is  an  indication  of  the  atmospheric  turbulence. 

The  capability  of  the  LDV  to  measure  turbulence  should  be  studied  in  more 
depth  and  further  exploited.  It  is  noted  that  the  LDV  can  measure  turbulence  at 
higher  frequencies  than  can  be  accomplished  by  anemometers.  Also,  the  LDV 
can  measure  velocity  fluctuations  over  space,  whereas  an  anemometer  can  only 
measure  velocity  fluctuations  over  time  at  a fixed  location. 

4.4.4  Double  VAD  Signature 

Figure  4-13  shows  a distinct  double  VAD  signature.  Two  distinct  signa- 
tures are  discernable,  and  it  appears  that  two  winds  are  being  measured.  The 
characteristic  is  not  anomalous  because  it  appears  on  many  of  the  signatures 
measured  at  approximately  the  same  time. 

A double  VAD  signature  has  appeared  previously  in  snow  (Ref.  6),  and  can 
be  plausibly  explained  in  snow.  For  the  present  condition,  a plausible  explana- 
tion is  not  clearly  evident.  The  phenomenon  could  be  wind  shear.  The  focal 
volume  covers  an  altitude  of  278  m jt25  m.  Wind  shear  may  easily  be  sufficiently 
large  to  cause  a difference  in  wind  speed  as  shown  in  Fig.  4-13,  especially  with 
the  small  quantity  of  turbulence  indicated  by  the  well-defined  sine  curves  in 
Fig.  4-13.  The  wind  shear  could  result  from  stratification  of  a very  stable 
atmosphere,  which  often  occurs  at  night.  Table  4-1  indicates  that  this  run 
was  taken  shortly  after  midnight. 
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FIGURE  4-13.  VELOCITY  AZIMUTH  DISPLAY  LINE-OF-SIGHT  VELOCITY 
SIGNATURE  SHOWING  A DOUBLE  SIGNATURE. 
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4.4.5  Data  Recorded  in  Rain 


On  23  June  1976,  data  were  recorded  in  rain.  In  the  analysis  of  data  in 
rain,  it  is  helpful  to  examine  the  characteristics  of  data  recorded  in  rain.  Figure 
4-14  shows  the  geometry  of  the  measurement  when  a significant  vertical  com- 
ponent is  imposed  on  the  airborne  particulates  (such  as  raindrops).  The  vertical 
component  is  shown  upward  to  coincide  with  previously  accepted  notation  for  the 
development  of  the  equations.  In  rain,  the  vertical  component,  w,  is  negative. 

For  laser  reflection  from  raindrops,  the  magnitude  of  the  line-of-sight  velocity 
component  is 

v^  = cos (t/ 2 + tan  ^ w/ v - a)  , 

where  \v^  + w^  is  the  speed  of  the  raindrops  (assuming  that  the  horizontal 
component  of  drop  velocity  equals  the  horizontal  wind  component),  and 

tt/  Z + tan  1 w/ v - a 

is  the  angle  between  the  velocity  vector  of  the  raindrops  and  the  line-of-sight 
direction  of  the  LDV.  For  small,  negative  values  of  w,  the  effect  of  the  vertical 
component  is  to  increase  the  line-of-sight  velocity  on  the  upwind  side  of  the  LDV 
and  to  decrease  the  line-of-sight  velocity  on  the  downwind  side  of  the  LDV.  On 
the  downwind  side,  the  velocity  component  of  the  raindrops  may  become  almost 
perpendicular  to  the  line-of-sight,  and  the  line-of-sight  velocity  component 
vanishes  on  the  downwind  side. 

The  beginning  of  the  rain  is  shown  sequentially  in  Figs.  4-15,  4-16,  and 
4-17.  Figure  4-15  presents  data  acquired  before  the  rain  began  and  shows  a 
normal  VAD  signature.  Figure  4-16  shows  both  reflections  from  natural  aerosol 
and  from  raindrops.  The  sign  convention  ol  the  derectified  line-of-sight  velocity 
signature  is  that  the  positive  portion  of  the  signature  is  on  the  upwind  side  of  the 
VAD  scan.  Therefore,  according  to  the  previous  paragraph,  the  highi  st  signa- 
ture is  the  signature  from  rain.  Figure  4-17  shows  data  in  ram.  Only  one 
negative  velocity  point  occurs.  Some  of  the  plots  have  no  negative  velocity 
points.  This  could  easily  result  from  the  trajectory  of  the  raindrops  being 
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approximately  perpendicular  to  the  line -of- s ight  direction  of  the  LDV  on  the 
downwind  side. 

It  is  postulated  that  the  cause  of  the  missing  data  was  that  the  frequency 
threshold  setting  was  set  too  high.  The  log  sheet  shows  the  setting  as  500  kHz 
(~2. 7 m/sec);  however,  only  two  points  below  4 m/sec  appear  in  Figs.  4-15, 
4-16,  and  4-17.  In  future  measurements  in  rain,  it  will  be  necessary  for  the 
operators  to  be  aware  that  rain  causes  a decrease  in  the  line-of-sight  velocity 
on  the  downwind  side  of  the  scan,  and  threshold  settings  must  be  made  accord- 
ingly. 


Three  additional  comments  are  relevant  to  wind  measurement  in  rain. 
First,  the  SEL-810  data  logging  program  used  for  the  rain  data  was  different 
from  that  used  in  the  rest  of  the  test.  After  the  test  described  herein,  it  was 
determined  that  the  data  logging  program  occasionally  skips  data  (about  one- 
half  second  skipped  for  every  third  scan).  The  problem  has  been  corrected, 
but  it  may  account  for  some  of  the  lack  of  data  for  negative  line-of-sight  veloc- 
ities. Second,  all  rain  data  were  recorded  at  the  444-m  altitude.  It  is  expected 

that  the  signal  strength  from  raindrops  may  be  stronger  at  a lower  altitude.  The 

4 

focal  volume  of  the  LDV  is  proportional  to  R (R  is  range).  The  laser  intensity 

-4 

at  the  detector  is  proportional  to  R .If  the  number  density  of  particles  (num- 
ber of  particles  per  unit  volume  of  atmosphere)  is  large,  the  number  of  par- 
ticles from  which  a scattered  signal  is  received  is  proportional  to  the  focal 

4 -4 

volume,  the  R term  cancels  the  R term,  and  the  returned  laser  intensity  is 
constant  with  changing  range.  This  is  the  condition  for  natural  atmospheric 
aerosol.  For  a low  number  density  of  particles  (as  for  rain),  the  number  of 
particles  in  the  focal  volume  at  a given  instant  in  time  is  small.  Increasing 
the  focal  volume  size  does  not  mean  that  the  number  of  particles  in  the  focal 
volume  will  increase.  Therefore,  the  returned  laser  intensity  may  decrease 
with  increases  in  range  for  a small  number  of  particles.  The  effect  of  range 
on  returned  laser  intensity  in  rain  is  a subject  for  further  investigation. 

Third,  it  is  known  that  rain  removes  aerosol  from  the  air.  This  causes 
no  problem  if  there  are  raindrops  from  which  the  laser  beam  can  be  reflected. 
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FIGURE  4-15.  VELOCITY  AZIMUTH  DISPLAY  SIGNATURE 
BEFORE  BEGINNING  OF  RAIN. 
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FIGURE  4-16.  VELOCITY  AZIMUTH  DISPLAY  SIGNATURE 
AT  THE  BEGINNING  OF  RAIN. 
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FIGURE  4-17.  VELOCITY  AZIMUTH  DISPLAY  SIGNATURE  DURING  RAIN. 


However,  the  adequacy  of  the  natural  aerosol  after  the  rain  stops  has  not  been 
evaluated.  From  the  Table  Mountain  tests  (Ref.  6)  the  aerosol  remaining  in 
the  air  after  a snow  was  barely  adequate.  At  lower  altitudes  the  quantity  of 
natural  aerosol  is  much  greater  than  at  Table  Mountain.  Therefore,  the  lack 
of  natural  aerosol  after  a rain  is  not  expected  to  be  a problem,  but  verification 
is  necessary. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  Lockheed  laser  Doppler  velocimeter  appears  to  be  able  to  measure 
wind  accurately  to  an  altitude  of  450  m.  A discrepancy  of  0 to  3 m/  sec  was 
observed  between  the  laser-measured  wind  speed  and  the  anemometer- 
measured  wind  speed.  The  discrepancy  occurred  at  all  altitudes.  The  dis- 
crepancy does  not  indicate  a deficiency  in  the  LDV  because  the  discrepancy 
occurred  at  altitudes  at  which  the  accuracy  of  the  LDV  had  previously  been 
verified,  and  a discrepancy  between  anemometer-measured  winds  and  winds 
indicated  by  double  theodolite  balloons  has  previously  been  reported  at  the 
WKY-TV  tower.  The  wind  speed  discrepancy  may  be  directionally  dependent. 

A conclusion  could  not  be  reached  because  almost  all  data  were  from  the  same 
direction.  Directional  agreement  between  laser-measured  winds  and  anemometer- 
measured  winds  was  good  with  almost  all  data  in  agreement  within  +10  deg  for 
3 -min  averages  and  within  +5  deg  for  15-min  averages.  Occasionally,  differ- 
ences of  +20  deg  were  observed. 

The  measurement  of  wind  in  rain  was  inconclusive.  Deficiencies  in  the 
data  were  noted,  but  very  simple  remedies  exist.  It  is  expected  that  a better 
signal  would  be  obtained  at  a lower  altitude.  For  wind  shear  measurement 
for  aircraft  approaches,  the  altitudes  of  interest  are  less  than  444  m.  The 
glideslope  altitude  is  approximately  444  m at  the  outer  marker,  and  it  is 
usually  the  region  inbound  from  the  outer  marker  which  is  of  most  concern 
for  wind  shear  related  to  aircraft  approaches.  A more  complete  measure- 
ment and  analysis  program  for  wind  measurement  in  rain  is  appropriate. 
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Appendix  A 

TYPICAL  VELOCITY  AZIMUTH 
DISPLAY  SIGNATURES 

Typical  VAD  signatures  are  presented  in  sets  of  three.  The  sets  may 
be  identified  by  the  run  number  and  time  given  in  the  legend.  The  first  figure 
of  each  set  is  the  raw  data.  The  second  figure  shows  data  which  have  been 
edited  and  filtered,  and  the  third  figure  of  each  set  shows  the  derectified  signa- 
ture (edited,  but  unfiltered). 
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FIGURE  A - 1 . TYPICAL  LASER  DOPPLER  VELOCIMETER  SIGNATURES 


ALTITUDE  It  KJB.O  lCTO» 

TflC  IS  l!'?»*57  OKLAHOMA  MOROf  RUN  ",  VAO  %/  t/TB  N0W1AN  *>175. 


?l  ROIKT  AVERAGE 


HORIZONTAL  ANOLX  IN  OCORETS 


FIGURE  A - 1 (Continued) 
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FIGURE  A- 1 (Continued) 
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FIGURE  A - 1 (Continued) 
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FIGURE  A-l  (Continued) 


A - 7 


ALTITUDE  IS  39.0  rCTtWS 

Tire  IS  IIM9*?9  OKLAHOMA  NOR03  VAO  6 / 9/76  NORMAN  HD1T5. 


HORIZONTAL  AHCAt  IN  OTOAICS 


FIGURE  A - 1 (Continued) 
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FIGURE  A-l  (Continued) 


A - 10 

1 


MOPfTONTM.  ANOLC  IN  DCOfltCS 


FIGURE  A - 1 (Continued) 
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FIGURE  A- 1 (Continued) 


A - 1 5 


I 


NMimra.  mi  m 


FIGURE  A 


1 (Continued) 


A - 17 


HORIZONTAL  AMOLC  IN  OCONOI 

FIGURE  A - 1 (Continued) 


HORIZONTAL  AMOLl  IN  OCOMZ* 

FIGURE  A-l  (Continued) 


ALTinjOC  It  ttt.o 

Tift:  It  *0,«*,I7 


k NOROS  NUN  s 


VAO  t/ll/lt  NCVtlAN 

et  POINT  AVCRAOC 


HORIZONTAL  ANOLC  IN  OUflUt 


FIGURE  A-l  (Continued) 
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FIGURE  A-l  (Continued) 
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Appendix  B 

TYPICAL  TABULATED  DATA 

Typical  tabulated  data  sheets  are  presented.  The  tabulated  data  show 
typical  comparisons  of  the  results  of  the  various  algorithms,  typical  values 
of  vertical  component  of  wind,  typical  values  of  the  second  and  third  harmonic 
for  the  Fourier  (spectral)  algorithm,  and  typical  values  of  the  standard  devia- 
tions of  the  various  parameters.  A complete  set  of  tabulated  data  exists  at 
the  Department  of  Transportation,  Transportation  Systems  Center,  Cambridge 
Ma  s sachu  setts . 
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TYPICAL  TABULATED  WIND  DATA  FROM  LASER  DOPPLER  VELOCIMETER 
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Appendix  C 

COMPARISON  OF  LASER-MEASURED  WIND  AND  TOWER-MEASURED 
WIND  FOR  1-,  3-,  6-,  9-,  12-  AND  15- MINUTE  AVERAGING  PERIODS 


FIGURE  C-l.  COMPARISON  OF  LASER-MEASURED  WINDS  WITH 
ANEMOMETER -MEASURED  WINDS. 
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Appendix  D 

COMPARISON  OF  LASER -MEASURED  WIND 
AND  ANEMOMETER -MEASURED  WIND  FOR 
3-  AND  15-MINUTE  AVERAGING  PERIODS 


Figures  for  the  15-min  averaging  period  are  eliminated  whenever 
the  run  is  sufficiently  short  that  there  are  less  than  three  15-min 
averages . 
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FIGURE  D-l.  COMPARISON  OF  LASER  -MEASURED  WINDS 
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Appendix  E 
REPORT  OF  INVENTIONS 


The  purpose  of  the  work  performed  under  tills  contract 
and  reported  herein  was  the  further  verification  of  the 
capability  of  an  existing  remote  sensing  device  to  measure 
atmospheric  winds  accurately.  It  was  desired  to  provide 
verification  of  the  accuracy  of  the  laser  Doppler  velocimeter 
to  higher  altitudes  than  the  highest  altitude  of  previous  tests. 

Because  the  purpose  of  the  test  was  the  use  of  established 
techniques  on  an  existing  device,  no  innovation,  discovery, 
improvement,  or  invention  was  made. 
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